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ABST RACT 

Calculations are presented for the 1999 opposition class mission and 
a procedure for obtaining similar occultation data for any other given 
Mars mission is given. Occultation data for a Mars orbiter in a 24.5 
hour parking orbit and a Mars base have been calculated for: sunlight 

occultation - the time in darkness; and radio communication occultation - 
the communication losses between the lander and the orbiter, the lander 
and Earth, and orbiter and Earth. 


CALCUL ATIONS 

Mars Orbiter Sunlight Occultation 

To find the time in darkness for a Mars Orbiting Spacecraft it is 
necessary to determine the orientation of the parking orbit with respect 
to the Sun. This is done by finding the angle between the semi-major 
axis of the orbit and the Mars to Sun line. This angle, a , (see 
Figure 1), is found using the following equation: 

a x = RAP - VE - P - L 


The values of RAP, VE, P and L are found using the trajectory data 
(1) and the Planetary Handbook (2). 

Once ot is known it is possible to find the points on the orbit 
corresponding to the beginning and end of occultation, defined as ^, and 

^ . This is done using the equation of the parking orbit and a 
2 

transformation between a reference frame centered on the orbit ( x, y) 
and one that lies along the Mars-Sun line (x^ , y ). These reference 
frames are shown in Figure 2. The second frame, (x^ , y^ ) is defined as 
a rotation of the orbit centered frame (x, y) through an angle ot 
followed by a translation d, where d is defined as: 

d - ae sin ot 

The quantity ae represents the distance from the center of the orbit 
to the focus of the orbit which is the center of Mars. 

The values of ^ and are found by finding the values of x 
that correspond to y^ equal to + the radius of Mars (RM), then con- 
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FIGURE 1 

GEOMETRY FOR DETERMINING MARSORBITER OCCU 
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FIGURE 2 

GEOMETRY FOR CALCULATING MARSORBITER OCCULTATIONS 




verting to the x,y frame. [Once ^ an H tit , , 

1 V 1 and V 2 are known,] Kepler's time 

equation then can be used to find the tl„e fro. perigee to * and * 

The difference in the two ti.es represents the duration of the 
occultation. 

The transformation between the two coordinate frames are found in 
the following manner: 

First, a rotation of the x,y frame through an angle a to yield the 
K l >y j axes (see Figure 1). 

X 2 = x cos a+ y S i n a 
* ~x sin ot + y cos^t 

Next, a translation along the y axis a distance d to yield the x", 
y axes. 

x = = x cos a + y sin a 

y yj - d = -x sin a + y cos ot _ ^ 

„ In . matrlX notatlon the transformation between the x,y frame and the 
X , y frame can be written as: 
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To find the values of x corresponding to y" 
the orbit must be found in terms of x ' and y " 
orbit in the x,y frame is: 


i RM the equation of 
The equation of the 


x 

a 


2 

2 



= 1 


or: 


xV ♦ y 2 a 2 - a 2 b 2 = 0 

using the transformation: 

x = x cos a - (y + d) sin ot 
y = x sin Qt + (y + d) cos ot 
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The equation of the orbit in terms 

• y O O O O 

x I b COS a + a sin u J \ 
+ [(y d ) 2 (b 2 sin 2 


of x and y" is: 

X " [2(y "+ d)(a 2 - b 2 ) 

O 0 O O 

a + a cos a) - a b ] = 


sin r *cos fi ] 


0 


Given y > = +RM the corresponding values of can be found using 
the quadratic formula, then these values of x^ and y" can be converted 
to the x,y frame using the transformation matrices. The values of ^and 
are found using the x and y coordinates of the orbit at the beginning & 
end of occultation. 


The 
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duration of 
to calculate 


beginning 

the occultation is found by 
the 1 ime from perigee at i/j 
3 


P 



( ty- e sin iff ) 


i* 


end 


using Kepler's time 
and ^ 2 . 


The duration of the occultation is: 


period of At once during 
the staytime as the longi- 
To obtain the minimum and 
on Mars arrival and depar- 


would be subjected to is 
highly dependant on the latitude of the landing site and the heliocen¬ 
tric longitude of Mars. The amount of daylight varies on Mars just as it 
does on Earth, since its equator is inclined to its orbit by 23.984 
degrees. To calculate the time in darkness, the following calculations 
are required. 

The geometry shown in Figures 3a and 3b represents the orientation 
of the axis of rotation of Mars with respect to the Sun. H represents 
the angle between the polar axis and the vertical as viewed perpendicular 
to the Mars-Sun line. Us uig Figure 3a, f) can be found: 


The orbiting spacecraft is occulted for a 
each orbit. The value of A t changes during 
tude of Mars changes and the orbit precesses. 
maximum occultations, At should be calculated 
ture. 

Mars Lander Sunl_i_ght Occ ultation 

The time in darkness that a Mars lander 
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3218—85 FIGURE 3 

GEOMETRY FOR FINDING 1 ? 
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CROSS SECTION A-, 







i = 23.984 DEG 
r - RM sin i 

y ■ M cos L = RM sin 1 cos L 
z = RM cos i 


tan r| = 


cos i 
cos L 


Figure 3b shows the orientation of Mars with respect to the Sun and 
a cross section of a particular latitude shows how Much of that latitude 
is in the sunlight. The time in darkness is found using the following 
calculations and Figure 3b. 


tan n 
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RM sin SL 
tan n 
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Communication Occultation Between The Mars Orbiter And Earth 
This calculation is performed using the same procedure as the Mars 
orbiter sunlight occultation except instead of a being used in the 
calculations, a different angle, 6 is used. 6 is defined as the angle 
between the semi-major axis of the parking orbit and the Mars to Earth 
line. The following calculations are required (see Fig. 1). 


Q - 


r 

e 


L - LE | 

a (1-e 2 ) ' 

1 + e cos© 


Barth 
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r 

■ 


Mars 



= r m ~ r e cos Q 

Y = r g sin Q 

r 

J = tan Y 
~X 

3 = a- J 


Once 3 is found the procedures for finding the communication occultations 
are identical to the sunlight occulations. Starting with equation 1, Bis 
substituted for a . 


Communication Occultatlon Between The Mars Lander and Earth 
The communications occultatlon between the Mars lander and the Earth 
is dependent on the same parameters that influence the sunlight occupa¬ 
tion. Since the plane of the Earth's orbit is inclined only 1.849 
degrees to the Mars orbit plane, the value of the duration of the 
communications occultatlon would be essentially equal to the duration of 
the sunlight occultatlon for a given latitude. 

Communication Occultatlon Between The Mars Lander and The Mars 
Orbiter 

The communication occultatlon of the Mars lander and the Mars orbi¬ 
ter is obtained by finding when the angle between the local vertical at 
the landing site and the position of the orbiter is greater than 90 
degrees. The geometry for this calculation is shown in Figure 4. It was 
assumed that the orbiter is directly above the lander when it is at 
perigee. The angle c represents the angle between the vertical and the 
orbiter. Communication is occulted as long as c is greater than 90 
degrees. Using Figure 4, c was obtained with the following calculations. 
Given a, e, y and ip. 

0-2 

t = tp 



2 IT RAD 


w = 


149 


24.5 


HRS 



3217-85 
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ORBITER 





$ = 2 IT — (D t 


J 

h 

k 


NOTE: $ is equal to this only when the parking 

orbit is retrograde, i.e., direction is opposite 
Mars's rotation. 

= | 0 - $ I 


= A7 + 

-i 

= sin 


2 

r - 2 RM r 
r sin J ~I 


cos j 



C = tt - K; h>r 


c = K ; h < r 

DATA FOR BASELINE. 1999 MISSION 

24.5 hour parking orbit 

r p = 3900 k “ = 36829.2 km 

a = 20364.63 kn 
e = .8084915 


MARS ORBITER SUNLIGHT OCCULTATION 



MARS ARRIVAL 

MARS DEPARTURE 

RAF = 

273.816 deg 

275.242 deg 

VE = 

-67.01 deg 

-67.01 deg 

P = 

335.323 deg 

335.323 deg 

L - 

168.915 deg 

195.73 deg 

a = 

16.588 deg 

8.801 deg 

At a 

.267 hr (16.06 

■in.) .264 hr (15 


MARS LANDER SUNLIGHT OCCULTATION 

No landing site has been baselined sc^ assume occultation of 1/2 
Mars day. 


At = 12.25 hrs 


COMMUNICATION OCCULTATIONS 
ORBITER TO EARTH 


L 

LE 

Q 


MARS ARRIVAL 
“ 168.91 

= 115. 

• 53.91 


MARS DEPARTURE 
195.73 
176. 


12.747 
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19.73 

73.477 



220.41 


193.588 


r 

e 

— 

.98367 

AU 

.99498 

AU 

r 

m 

= 

1.6604 

AU 

1.62526 

AU 

X 

c 

1.0809 

AU 

.68449 

AU 

Y 

= 

1.315 

AU 

.33589 

AU 

J 

= 

50.59 deg 


26.138 deg 


3 

= 

34.002 deg 


17.337 deg 


At 

= 

.2823 hr 

(16.93 min.) 

.2677 hr (16 

. 06 min.) 

LANDER TO EARTH 




Since 

no 

landing site 

has been chosen, 

assume communication occul 


tatlon occurs for half of the Mars day. 


AT = 12.25 hr 

LANDER TO ORBITER 

These calculations were made by calculating c for = o to 360 
degrees. The results are shown below: 


tp (hr) 

^(deg) 

c (deg) 


0 

0 

0 



.3 

18 

79.2 



.4 

24 

108.7 

occultation 

begins 

7.5 

140 

89.2 

occultation 

ends 

12.25 

180 

0 



17.0 

220 

89.2 



CM 

t- 

r-f 

222 

93.2 

occultation 

begins 

24.1 

336 

108.7 

occultation 

ends 

At = 7.1 hrs 

occuring twice 

every orbit. 
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LIST OF SYMBOLS 


a = semi-major axis of an orbit 
b * semi-minor axis or an orbit 

c = angle between local vertical at landing sight and the Mars orbiter 
d = distance from origin of x,y reference frame to the origin of the 
x 1 f frame 
e * orbit eccentricity 

i = inclination of Mars’ equator to the ecliptic 

1 = latitude of landing site 

L = heliocentric longitude of Mars 

LE = heliocentric longitude of Earth 

P = heliocentric longitude of Mars perihelion 

RM = radius of Mars 

RAP = right ascension of perigee of Mars parking orbit 

r = distance from center of Mars to orbiter 

r = distance from center of sun to center of Earth 
e 

r = distance from center of sun to center of Mars 
m 

tp 1 = time from perigee that occultation begins 

tpg = time from perigee that occultation ends 

A t = duration of occultation 

VE = heliocentric longitude of Mars vernal equinox 
x,y = coordinates in the x,y reference frame (Mars parking orbit) 
x* , y' - coordinates in the x'* , y " reference frame 
a = angle between semi-major axis of parking orbit and Mars-Sun line 

0 «= angle between semi-major axis of parking orbit and Mars-Earth line 

ip = eccentric anomaly 

n = angle between vertical and the Mars polar axis as seen perpendicular 
to the Mars-Sun line 
p = gravitational constant kmg/seCg 
0 * true anomaly 
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